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Abstract 

Following the original idea of Debye, we define and extract a gauge-invariant 
screening mass from the complex static in-medium heavy-quark potential Vqq, 
recently obtained from lattice QCD. To this end we derive a field theoretically 
motivated analytic formula that faithfully reproduces both the screened real¬ 
as well as the imaginary part of the lattice potential with a single temperature 
dependent fit parameter mu{T). Using values of the real part of Vqq in a 
gluonic medium, we obtain Debye masses compatible with predictions from 
HTL perturbation theory. 
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1. Introduction 

The concept of a screening mass helps us to intuitively understand the in¬ 
tricate interplay of interactions that take place, once a test particle is inserted 
into a medium of charge carriers. Originally Debye and Hiickel [1] investigated 
the behavior of ions in electrolyte solutions and found that their interactions 
could be understood by an exponential suppression of the vacuum Coulomb 
potential due to the presence of classical Boltzmann distributed charges. This 
in-medium modification, translated into the modern language of quantum field 
theory, amounts to a medium induced dressing of the mediating gauge bosons, 
bestowing them with an otherwise absent mass. In the context of perturbative 
quantum electro-dynamics this mechanism of a thermal mass generation is well 
understood, as arising from the zero momentum limit of the gauge invariant 
longitudinal gluon self energies. In perturbative QCD, only the leading order 
and the logarithm at next to leading order (NLO) of the Debye mass can be 
computed, the constant at NLO already receives non-perturbative contributions 
[5]. On the lattice, its definition encounters another difficulty, unlike in QED 
where the Debye mass can be extracted from the electric field correlator, in 
QCD the electric field itself is not gauge invariant. Several approaches were 
proposed to circumvent this problem using e.g. effective theories obtained by 
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dimensional reduction [3] , spatial correlation functions of gauge invariant meson 
correlators |3] or the behavior of the color singlet free energies Enam. 

Here we return to the original idea of Debye and identify a physical observ¬ 
able, the static heavy-quark potential V (r) between a color singlet configuration 
of a quark and anti-quark, to non-perturbatively define a screening mass for 
QCD. In the vacuum the potential exhibits both a perturbative Coulombic, as 
well as a non-perturbative string-like behavior jS]. Both of these features will 
receive in-medium modification, as can e.g. be seen in a recent lattice QCD 
determination [9]. An additional complication arises from the fact that the po¬ 
tential is in general a complex quantity m. due to the presence of scattering of 
light medium degrees of freedom with the color string spanning in between the 
heavy quark and anti-quark. A meaningful description of the relevant physics 
must therefore necessarily capture both the effects of screening of the real part 
of the potential ReH(r) as well as Landau-damping related to ImH(r). Our 
strategy hence is to find a field-theory motivated parametrization of the poten¬ 
tial that depends only on a single temperature dependent parameter TO£)(T), 
which we will be able to identify with the Debye screening mass. 

In the literature two contributions along this path can be found. On the 
one hand, Ref. [7] proposed an analytic function for the real part of a medium- 
modified Cornell-type potential by working within the fully classical setup of 
a test charge surrounded by Boltzmann distributed charge carriers, introduced 
originally by Debye and Hiickel [T] and generalized to non Coulombic potentials 
in Ref. m- While it was shown that the resulting parametrization of Re[V] 
can reproduce the lattice data quite well, it required the introduction of a sec¬ 
ond temperature dependent fit parameter which complicates the interpretation 
of the screening in terms of a unique screening mass. Since this classical ap¬ 
proximation is unable to accommodate an imaginary part of the potential, no 
statement about Im[V] was made. On the other hand an interesting approach 
was proposed in Ref. [12] , which attempts to describe the screening of the poten¬ 
tial in terms of actual in-medium field theory. The authors make the assumption 
that the in-medium potential arises from the vacuum potential by multiplying it 
with a field-theory determined complex permittivity in momentum space. Using 
the perturbative permittivity calculated in the HTL approximation it was pos¬ 
sible to reproduce the known real- and imaginary part of the potential in HTL. 
Unfortunately applying the permittivity to the linearly rising potential lead to 
unphysical results. The real-part does not decay exponentially but retains an 
^ l/r behavior, which does not describe the lattice potential and hints at the 
fact that the screening of the vacuum potential is not captured self-consistently. 
Even more severely, the resulting imaginary part, which from physical considera¬ 
tions must go to a constant at large distance, namely twice the Landau damping 
of a single quark, diverges logarithmically. 

Our study combines the strength of these two approaches bringing together 
the concept of a generalized Gauss law from Ref. m with the characterization 
of in-medium effects through the perturbative HTL permittivity. The use of 
the Gauss law, a non-local concept, leads to a self-consistent description of both 
screening and damping effects evading the unphysical behavior of Ref. [12]. 
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2. The static in-medium inter-quark potential 


The static potential acting between a color singlet pair of a quark and anti¬ 
quark immersed in a thermal medium of light quarks and gluons constitutes the 
basis for our gauge invariant QCD screening mass. The definition of the heavy 
quark potential at finite temperature itself has been a long standing problem 
in thermal field theory. Using an effective field theory description m, based 
on the separation of scales between the mass of the heavy quark and its typical 
momenta, it had been worked out how to define the potential in vacuum from a 
dynamical QCD observable, the real-time Wilson loop W{t,r). With the matu¬ 
ration of modern effective field theories [Hj, such as NRQCD and pNRQCD it 
became possible to take into account the presence of the additional scale repre¬ 
sented by temperature and to systematically extend the validity of the potential 
definition to finite temperature |15j . It reads 


V (r) = lim 

t—^OO 


idtW{t,r) 

W{t,r) 


( 1 ) 


This potential was first calculated in hard thermal loop (HTL) resummed per¬ 
turbation theory |10j and was found to be complex valued: 


Thtl(»’) = — ds 


^-mor 

mo -\ -h iT(j){mor) 


+ Oig^), 


( 2 ) 


where 

( 3 ) 

and a factor Cp has been absorbed in the definition of the coupling constant 
ds = 4 ^^ to match the literature on phenomenology. Its real-part indeed 

shows the typical Debye screened behavior. A non-perturbative, i.e. lattice 
QCD based, determination of the potential however remained a conceptual and 
technical challenge, which has only recently been overcome in a satisfactory 
way [i m El na [n]. The central hurdle is related to the fact that lattice 
simulations are performed in Euclidean time and have no direct access to real 
time quantities such as W (r, t). 

A possible way around this limitation was proposed in Ref. |20j with a first 
attempt at an implementation presented in Ref. m- The underlying idea is to 
use a spectral decomposition of the Euclidean Wilson loop W (r, r) to relate the 
Euclidean an Minkowski time domain 


W{T,r) 


j due o 


J due 


W{t, r). 


This equations can be combined with Eq. Q to define the potential in terms of 
the Wilson loop spectral function 


V{r) = lim 

>-oo 


duue '‘^*'p(u,r)l 


due '-‘^*p{u,r). 


(4) 
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At that stage the definition requires precise knowledge of the spectrum of the 
Euclidean Wilson p(uj,r), which can in principle be obtained from an inverse 
Laplace transform of datapoints VE(r„,r), n = l..Nr simulated in lattice QCD. 
Bayesian inference plays an important role to give meaning to this otherwise 
ill-posed problem. In practice however carrying out the inverse Laplace trans¬ 
form posed a formidable challenge to standard methods, such as the Maximum 
Entropy Method m or extended MEM [22]. In fact, in the end it required the 
development of a novel Bayesian inference method [19j . 

Even if the Wilson loop spectrum can be determined reliably from the lattice, 
a second difficulty lies in the infinite time limit in Eq. i)- Using only the 
symmetries of the real-time Wilson loop, it was shown in Ref. m that the 
physics of the potential manifests itself in the lowest lying peak in the spectrum, 
which takes the shape of a skewed Lorentzian 

|ImU(r)|cos[Recroo(?')] — (ReU(r) — a;)sin[Re(Too(c)] 

ImU(r)2-p (ReU(r) — a;)2 

-I- Coir) + ci(r)(ReU(r) - w) -I- C 2 (r)(ReU(r) - .... 

Inserting this functional form of the spectrum into the definition it was 
conhrmed that the position and the width of this peak encode the values of the 
real- and imaginary part of the potential [23] . 

To put this extraction strategy into practice we generated anisotropic = 
3.5 quenched QCD configurations at /? = 7 with temporal extend = 24 ... 96, 
i.e. spanning T = 839... 210MeV (Tc ~ 271MeV), and extracted both the real 
and imaginary part of the potential as reported on in Ref. [Sj. Eor the current 
work we added an additional set of A^conf = 900 low temperature configurations 
at Nr = 192, i.e. T = lOSMeV, that will be used in the determination of the 
Debye mass. 

3. An analytic parametrization of the heavy-quark potential 

The starting point of our derivation is the generalized Gauss law introduced 
in Ref. m for an electric held the form E = qr°- 

^ ( 5 ) 

Using the relation —VU(r) = E(r) Eq. [^reduces to the well known expression 
for the Coulomb potential for a = —l,q = Os, [og] = !> while a linearly rising 
potential corresponds to a = 1 ,9 = cr, [cr] = GeV^. 

Let us Hrst have a look at the original argument by Debye and Hiickel. In 
their work the above equation is modihed through the presence of a background 
charge density (p(r)) 

^ = 47r(^(rQ + (p(rO)), (6) 
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which represents a Boltzmann’s distribution for the charge carriers at temper¬ 
ature T = 1//3, 

(p(^) = Q - no) - g - no) , (7) 

where in our context the hrst term stands for particles and the second for an¬ 
tiparticles. no is the charge density in the absence of the test charge. If the 
resulting in-medium potential is weak, we can expand the exponential in Eq. 
as 

{p{r^) =-2qPnoV{f). ( 8 ) 

When plugged into Eq. Q, we obtain HU: 

= 47r9<5(r0, (9) 

where A = SirquoP. The appearance of the term with prefactor A is a manifes¬ 
tation of the linear-response character of this approximation. For the Coulombic 
part of the potential {a = — l,q = ds) we have 

-VVc(r) -h AcVc = 47rd,(5(70. (10) 

On the other hand for the string case {a = \,q = a) one finds 

-^—fi^+AsVs{r)=^T:<T5{r), (11) 

where As = STrcrno/S. If no is indeed the unmodified charge density it has to 
take the same vale for both cases, which will allow us to relate Ac and As in 
the following. 

Now let us return to the generalized Gauss law of Eq. [^for a Coulomb charge 
(a = — 1 ) written in momentum space 

P^Vcip) = dTrds. (12) 


In the language of held theory, the effects of the medium is to dress the point 
charge with a cloud of charge carries. It amounts to washing out the coordinate 
space delta function on the RHS with an in-medium permittivity e(p, rn/j), 
which encodes the full temperature dependence of the problem: 


P^Vcip) = 47r- 


0-8 


e{p,mD)' 

Similar to Ref. [12], we use the perturbative HTL expression 


-ir^ ^ ■ rj. Pmjj 


(13) 


(14) 
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henceforth assuming a medium of weakly coupled quarks and gluons, in which 


our test charge is immersed. Inserting this formula into equation (13) and 
multiplying by ^ , we obtain: 


P^Vcij^ + rnl)Vc{p) = 47ras(l - it^T 


'-n 


p(p2 +m%) 


(15) 


The inverse Fourier transform of the real part of Eq. exactly reproduces 
the linear-response expression of (10), allowing us to identify Ac = rnjj and in 
turn gives an expression for the charge density 


no = 


od-‘ 


Sttq;, 


(16) 


The imaginary part arising on the RHS of Eq. (15) can also be Fourier trans¬ 


formed to coordinate space, which completes our generalized formula for the 
in-medium modification of a Coulombic test charge 


-V‘^Vc{r) + m%Vc{r) = ag ^47ri5(r) - iTmjjg{mDr)'j , (17) 


with 


g{x) = 2 f 
Jo 


dp 


sin(pa:) p 


(18) 


/o p^ + l 

The solution to Eq. 0 with the physical boundary condition R.eFc'(r)|^^^ = 
0, ImVc'(?')lr.=o = 0 '5rImFc(’')lr=oo = coincides with the HTL potential 
obtained by Laine et al. uni 


-l-mDr), lmVc{r) =-asT(j){x), (19) 


ReV'c’(»’) = -- 


where 4’{x) has been defined in Eq. ([^. 

Now we turn to the string-like, i.e. linearly rising part of the test-charge 
potential. Note that in contrast to the Coulomb case, for a = 1 we cannot 
transform Eq.([^ into a simple relation for the Fourier space potential akin to 
Eq.([T^. Instead, as a first step, we return to Eq. ([TT|) to obtain an explicit 


expression for the coefficient Ag. Taking the charge density no to be the same, 
irrespective of the test charge being Coulombic or string-like one then obtains 


together with Eq. (16) 


Ag = pL = m 




Ctg 


( 20 ) 


In the context of the linear-response pitcture this tells us that the strength of the 
in-medium modification of the linear potential is controlled by the parameter p 
and not mo itself. To arrive at the defining equation for the medium modified 
string-potential, we assume the validity of the linear response approximation if 
the charge density on the RHS of Eq. ( [TT] ) is replaced by the medium modified 
charge density obtained from HTL in Eq. (lli 


^ ^ = a(^dTT5{r) - iTm\ig{mDr)^- ( 21 ) 
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Note that these expressions differ significantly from the ones used in [7], where 
As was chosen on purely dimensional grounds to be rn\). As we will see in the 
next section, our choice motivated by the charge density allows one to adequately 
fit the lattice data as is and thus makes the introduction of any other fitting 
parameter unnecessary. 


Solving for the real part of Eq. (211 with the same physical boundary con¬ 


ditions as in the Coulomb case we find 

r[i] 


ReVs{r) = + 

24 .y/TT M 


liil 

2r[|] 




( 22 ) 


which besides the definition of /i also differs in a factor y/2 in the argument of 
the parabolic cylinder function D^{x) compared to Ref. [7]. The imaginary part 
of the in-medium modified string potential can be written in a closed form as a 
Wronskian solution: 


ImVs(r) = —i 


am 


D 


T 




= -idsT'0(/rr), 


(23) 


with 


iIj[x) = D_ii2{V2x) / dyYleD_ii2{iV2y)y‘^g{ymD/ 

Jo 

poo 

+ReD_i/2{iV2x) / dy D_i/2{V2y)y'^g{ymD/n) 

J X 

pOO 

-^- 1 / 2 ( 0 ) / dyD_x/2(,'/2y)y^g[ymD/y), 

Jo 


(24) 


This completes our derivation of the analytic expression for the real- and imag¬ 
inary part of the in-medium QCD heavy-quark potential. 

Let us have a look at the behavior of the solution just obtained. In the limit 
of zero temperature, i.e. vanishing mo, we recover the Cornell potential in 
the real part. As expected, at small distances the Coulombic real part behaves 
as 1/r whereas the string shows a linear rise with r. The imaginary part on 
the other hand rises according to for the Coulombic part and with for 


the string. I.e. the Coulombic HTL part (191 dominates at small r. At large 


distances we again find that the Coulombic part dominates the real part and 
behaves just like the naive Debye screened potential exp(—m]jr)/r. The fact 
that the string part dies off much more rapidly as exp(—TO|,r^/2) is the reason 
why we can actually identify the parameter mo with the Debye mass, when 
fitted to the functional form of the lattice potential. At asymptotically large 
distances both the Coulombic and string imaginary part saturate to a constant 
as required. 


4. Determining the Debye Mass through the complex in-medium 
heavy-quark potential from lattice QCD 

Our goal is to use the derived analytic expression for the in-medium poten¬ 
tial to extract the Debye mass from the static inter-quark potential recently 
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measured in lattice QCD. In this work we focus on the case of a purely gluonic 
medium, for which both Rey, as well as ImV have been determined at various 
temperatures in the phenomenologically relevant region around the deconfine¬ 
ment transition. Since only a single fit parameter mjj is required, we will carry 
out fits solely to the real-part of the potential, so that the agreement or dis¬ 
agreement of the corresponding imaginary part can serve as a crosscheck of our 
approach. 

We assume the values of the strong coupling ds and string tension a not 
to vary with temperature T, as they characterize the properties of the test 
charge to be inserted in the medium. Their values hence have to be determined 
in vacuum. In the absence of a true T = 0 lattice measurement, we use the 
newly generated lattice ensemble at T = 105MeV instead and fit the small to 
intermediate r region of ReR. There lattice cutoff artifacts are expected to be 
small, the lattice determination of the real-part is most reliable and the effects 
of the small but finite temperature are not yet significant. The particular nature 
of the lattice normalization of the potential introduces a constant shift c, which 
we also determine 

5^ = 0.206 ±0.011, cr = 0.174 ±0.011GeV^ c = 2.60 ± 0.023GeV. (25) 

Varying the fitting range up to a maximum of six steps at the upper and lower 
end of the fitting interval yields the error estimates shown. The only remaining 
free parameter at finite temperature then is the Debye mass niD- As can be seen 
by the agreement of the solid lines and data points in the left panel of Fig{^ its 
tuning alone allows us to achieve an excellent fit of the real-part of the potential 
at all temperatures, both qualitatively and quantitatively. To account for the 
propagation of the error on the low temperature parameters, besides changing 
the fitting range on the finite T potential, we also use in each range different 
combinations of the values for Q:s,cr and c according to the uncertainties from 
the fits at T ~ 0. The values we obtain for the Debye mass together with their 
error estimates are given in Tabj^ 

Note that while we have determined mo solely from an inspection of the 
real-part, the resulting values for ImV also agree reasonably well with the lat¬ 
tice data (Figj^ right). At high temperatures and small distances, where the 
lattice reconstruction of the potential is most reliable, we even find quantitative 
agreement with the analytical form within statistical errors. At larger distances 
we expect that the lattice data-points are indeed larger than the actual values 
of ImV, as the underlying extraction from spectral widths leads to unphysically 
large values due to a diminishing signal to noise ratio. 


SU(3): Nr 

24 

32 

40 

48 

56 

64 

72 

80 

96 

192 

r[MeV] 

839 

629 

503 

419 

360 

315 

280 

252 

210 

105 

N 

' meas 

3270 

2030 

1940 

1110 

1410 

1520 

860 

1190 

1800 

900 

molMeV] 

852±60 

709±70 

654±62 

537±59 

444±52 

264±62 

309±86 

401±62 

328±78 

0 


Table 1: Quenched lattice parameters {0 = 7,^i, = 3.5, Ns = 32) and values of the extracted 
Debye mass 
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Figure 1: (left) One parameter fit (solid lines) to the real part of the potential (data points) 
obtained in quenched QCD. (right) Imaginary part of the potential in quenched QCD (data 
points) and the values obtained from our analytic expression using the Debye mass fitted in 


ReV. 


Close to Tc at T = 315MeV our analytic postdiction of ImV appears to 
lie rather far away from the lattice data. The reason for this, we believe, is 
twofold. On the one hand the extracted values for ImV at this temperature are 
rather imprecise, with the statistical errorbars still sizeable. On the other hand 
the determination of mu from ReF also seems to lack accuracy. Indeed if we 
inspect the corresponding fit curve on the left of Figj^ we can see already by 
eye that is too steep compared to those at neighboring temperatures. Using a 
higher statistics estimate for uid from ReV at temperatures around the phase 
transition in the future, will most certainly yield closer agreement between the 
analytic expression and the actual lattice imaginary part. 

5. Debye mass and its match to HTL 

The temperature dependence of the extracted Debye masses can be com¬ 
pared to predictions from resummed HTL perturbation theory. According to 
Ref. [2] the Debye mass at leading log order can be written as: 



+cTg{fiT)^ + dTg{fiT)^ 


(26) 


where /r denotes the renormalization scale, g the running QCD coupling, and c 
and d are constants that represent non-perturbative contributions which need to 
be determined from a fit to the data. As there exists a clear degeneracy between 
the variation of the renormalization scale and the variation of the higher order 
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Figure 2: (left) Debye masses (blue) and HTL fit (red) of the temperature dependence ob¬ 
tained on quenched lattices. Error bars includes both the statistical error, as well as those 
from changing the fitting range and propagation of the error of ds and a determined at low 
temperature, (right) The dimensionless ratio of for use in phenomenological model¬ 

ing. 


contributions, parametrized by c and d, we choose to fix /r according to the 
usual convention = ttT and fit c, d from the obtained Debye masses in the 
previous section. 

For the running of the coupling (/(/r) we utilize the four loop result of 
Ref. [23] setting Kqqd = 0.216, appropriate for quenched QCD. The fit yields 
c = —0.40 ± 0.06 and d = 0.21 ± 0.06, which is shown in Fig. Note that 
the values obtained for c and d are small, which implies quite good agreement 
between hard thermal loop perturbation theory and the lattice extraction, even 
at the low temperatures probed here. In fact, one could argue that for a HTL 
fit, only the high temperature points should be used, for which a perturbative 
expansion is expected to apply. We see however that including the lower tem¬ 
peratures close to Tc does not change the determination of the parameters c, d 
in an significant way. 

In previous studies in quenched QCD, Debye masses were e.g. obtained by 
fitting not the proper heavy-quark potential but the color-singlet free-energies 
with a simple Coulombic Debye-screened form |5] . Comparing, we find that the 
values obtained here lie consistently lower than these previous estimates. 

We would like to note that even though the lattices used to determine the 
heavy-quark potential deployed in this study are quite finely spaced, no contin¬ 
uum extrapolation has been carried out. Therefore it is doubtful whether the 
values of the Debye mass shown on the left of Figj^ together with the values for 
ds and tr can be directly used in phenomenological models in the continuum. 
As a possible remedy we follow Ref. [7| in providing the ratio of TO£)/-\/ct on 
the right of Figj^ in which some of the systematic uncertainties arising from a 
finite lattice spacing might be expected to cancel. 
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6. Conclusion 


Based on a combination of the generalized Gauss law, introduced in m and 
the in-medium modification of its point charge distribution by a weakly inter¬ 
acting medium of light quarks and gluons, described by the HTL permittivity, 
we derived an analytic expression for the real- and imaginary part of the static 
inter-quark potential at finite temperature. After fixing the strong coupling 
and string tension at low temperature, we are able to qualitatively and quan¬ 
titatively reproduce the real-part of the potential measured in lattice QCD by 
fitting a single temperature dependent parameter m£)(T), which is proposed as 
gauge invariant screening mass in QCD. The temperature dependence of mjj we 
obtained in a purely gluonic medium agrees well with the predictions of HTL 
perturbation theory. Using the fitted values for mo we furthermore find that 
a quite successful postdiction of ImU at high temperatures is possible. Agree¬ 
ment with ImU at smaller temperatures, currently hampered by uncertainties 
in the fit of to the real-part, should improve once higher statistics has been 
collected. 

We hope that phenomenological modeling will benefit from the derivation of 
a well motivated and lattice data validated parametrization of both ReU and 
ImU. In addition our study opens up the possibility to extract the imaginary 
part of the potential in full QCD simulations, in which up to date only the real 
part has been determined in a reliable fashion. 
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